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Abstract 



For the orthogonal- unitary and symplectic-unitary transitions in random matrix the- 
ory, the general parameter dependent distribution between two sets of eigenvalues with 
two different parameter values can be expressed as a quaternion determinant. For the 
parameter dependent Gaussian and Laguerre ensembles the matrix elements of the deter- 
minant are expressed in terms of corresponding skew-orthogonal polynomials, and their 
limiting value for infinite matrix dimension are computed in the vicinity of the soft and 
hard edges respectively. A connection formula relating the distributions at the hard and 
soft edge is obtained, and a universal asymptotic behaviour of the two point correlation 
is identified. 



1 Introduction 

Random matrix ensembles for quantum Hamiltonians are defined in terms of the con- 
straint imposed on the general Hermitian matrix representing the Hamiltonian by time 
reversal symmetry and spin- rotation invariance (when relevant). This leads to ten distinct 
random matrix ensembles [p]]. The corresponding eigenvalue probability density function 
(p.d.f.) can be computed exactly when the independent elements are chosen to have a 
Gaussian distribution. One then finds that the ten distinct random matrix ensembles 
belong to one of six different eigenvalue p.d.f. 's. These are 

1 £ a-2 



i = i i<j<k<N 
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which defines the Gaussian random matrix ensemble, or 

i N 

^rUe-^r 1 n 1*2-* *<>° (i-2) 

l=i l<j<k<N 

which defines the Laguerre random matrix ensemble, where (3 = 1, 2 or 4. 

The study of these random matrix ensembles can be generalized to include a parameter 
r. This allows the situation to be studied in which the system evolves between two distinct 
eigenvalue p.d.f.'s as r is varied. For the Gaussian ensemble this can be achieved by 
choosing the joint distribution of the elements of the random matrix to be 

P(X (0) ; X; t) = A p>T exp ( - (3Tt{(X - e~ T X (0) ) 2 }/2|l - e - 2r |). (1.3) 

In this expression X is a real symmetric {(3 = 1), Hermitian {(3 = 2) or self dual quaternion 
(P = 4) X x X matrix, and X® is a prescribed random matrix which must belong to a 
subspace of the symmetry class of X. For r — ► oo the p.d.f. is independent of X^ ) and 
is that of the standard Gaussian ensembles, with p.d.f. ( [I . 1|) . For general r Dyson 
proved that the eigenvalue p.d.f. p = p(xi, . . . , x^; r) satisfies the Fokker-Planck equation 

with 

1 N 

W = WW = -J2 X J~ E loglxfc-^-l (1.5) 

i=i i<i<fc<^ 

(the superscript (H) is used because of a connection with the Hermite polynomials), 
subject to the initial condition that p agrees with the eigenvalue p.d.f. of X^ ^ at r = 0. 

For the Laguerre ensemble, a parameter dependent theory can be developed (see 
e.g. 0]) by considering the non-negative matrices A = X^X, where X is an n x m 
parameter-dependent matrix 

P(XW; X; r) = A PtT exp ( - /3Tr{(X - e^X^^X - e - r X(°))}/2|l - e - 2r |) (1.6) 

(c.f. ( |1.3|) ) and the elements of X are real (/3 = 1), complex (/3 = 2) or quaternion real 
([3 = 4), and X^ -* is a prescribed random matrix which must belong to a subspace of the 
symmetry class of X. For r — >• oo the corresponding p.d.f. for the square of the eigenvalues 
of A is given by ( |1.2|) with a = n — m + 1 — 2//3, while for general r the eigenvalue p.d.f. 
p = p(a; 2 , . . . , x 2 N ; r) satisfies (|1.4j) with 

i N N 

W = W^ = -Y / x 2 J -(a + l/P)^ogx j - J2 Iog|^-^| (1.7) 
^ i=i j=i i<j<k<N 

(the superscript (L) is used because of a connection with the Laguerre polynomials). 

Associated with the parameter dependent eigenvalue p.d.f.'s are parameter dependent 
distribution functions 

P(n+m)(%l j • • • ) X n^'l T U x l \ ■ ■ ■ i X m ') r 2) (1-8) 
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The distribution function p( n+m ) is defined such that the ratio 

P(n+m){ x i i ■ ■ ■ i x n ^ j Tlj j • • • j x \ni T 2) / P(n+(m-l)){ x l i ■ ■ ■ i x h^'i r li x ^ \ ■ ■ ■ i ^m-li T "i) 

gives the density of eigenvalues at x$ when the parameter equals r 2 , given that there 
were eigenvalues at x± \ . . . , x^ when the parameter equalled r 1; and that there are also 
eigenvalues at x^\ . . . , x$ with parameter r 2 (it is assumed r 2 > t%). 

For the parameter dependent Gaussian and Laguerre ensembles, quaternion determi- 
nant expressions have been derived for p( n+m ) in the case that (3 = 2 in ( |1.4| ) and initial 
condition proportional to with (3' = 1 and 4 |7|] . Because each (3 corresponds to 

a different symmetry class of the Hamiltonian (orthogonal, unitary and symplectic sym- 
metry for f3 = 1,2 and 4 respectively), this describes the transition from orthogonal or 
symplectic symmetry to unitary symmetry as the parameter r is varied from to oo. 

The determinant expressions apply for the finite system. As revised in Section 2, they 
involve skew orthogonal polynomials known from the study of the distribution functions 
for r = in the (3 = 1 and 4 cases || ||, the orthogonal polynomials related to the 
r = 0, (3 = 2 case, as well as the elements of a connection matrix which specifies the 
orthogonal polynomials in terms of the appropriate skew orthogonal polynomials. In 
the Gaussian case, this matrix can be written down from knowledge of the expression 
for the skew orthogonal polynomials in terms of the orthogonal polynomials (the Hermite 
polynomials). In the Laguerre case this step requires some calculation which is undertaken 
below. 

Once all these quantities are specified, the general (n + m)-point distribution function 
in the finite system is known explicitly. The immediate task is to compute the scaled 
thermodynamic limit. This is done in Section 3 for the Laguerre ensemble in the neigh- 
bourhood of the hard edge (the smallest eigenvalues, which are constrained to be positive), 
and in Section 4 for the Gaussian ensemble in the neighbourhood of the soft edge (the 
largest eigenvalues). 

In Section 5 inter-relationships between the distribution functions are discussed. This 
is done by analyzing certain limits in ( |1.8| ): the limit T\,T2 — > oo, t\ — r 2 fixed; the 
a — ► oo limit of the hard edge; and the bulk limit corresponding to large distances 
from the edge. Also considered is the asymptotic form of the dynamical density-density 
function pT l+1 JX,Y;t), which is shown to exhibit a universal non-oscillatory decay for 
X, Y, t — > oo. A consequence of this decay is an 0(1) fluctuation formula for the variance 
of a slowly varying linear statistic. In the final subsections a realization of the density 
at the soft edge for (3 = 1 initial conditions is given by an empirical computation based 
on the parameter dependent matrices ( |1.3| ) with a real symmetric matrix, and the 
explicit form of the static distribution functions at the hard and soft edges for /3 = 1 and 
4 is noted. 

In the Appendix it is shown how our results for the Laguerre ensemble with iV finite 
coincide, in the static limit r — > 0, with formulas given recently by Widom M. 
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2 Revision of the formalism 



At the special coupling (3 = 2, the Fokker-Planck equation ( |1.4| ) with W given by ( |1 . 5| ) or 
(|1.7| ) can be transformed into an imaginary-time Schrodinger equation for N independent 
fermions (see e.g. H). As such the exact iV-body Green function G({x^}, {xj}; r) can 
be written as a determinant involving the corresponding single particle Green functions. 
In terms of the new variable t/j = Xj in (|1 . 1| ) and yj = in ( |1.2|) we have previously 
shown H that 

G(a;[ 0) , ■ ■ ■ ,xP;xi, - ■ ■ ,x N ; r)dxi ■ ■ ■ dx N = G(y[°> , ■ ■ ■ ,y^';y 1 , ■ ■ ■ ,y N ; r)d yi ■ ■ ■ dy N (2.1) 
where 

ni (°) (°) ^ 



N 

,E r/2 TT 



w (Vj) TT (Vk-Vj) , ,r / (0) \ 1 , on x 

— tk II -jo] — jph det ig(y) %yfe;-r)] ijfc= i,..,Ar (2.2) 



Here 



g(y^,y;r) = ^w(y^)w(y)J2pAy m )pAy>~ ljT (2-3) 

j=0 

and {pj(x)} are the orthonormal polynomials associated with the weight function w(x). 
In the Gaussian case 

w(y) = e-y\p n ( y ) = 2- n /\- 1 /\n\y 1 / 2 H n (y), ln = n+ 1 -, (2.4) 

where H n (y) denotes the Hermite polynomial, while for the Laguerre ensemble 

r(n + l) u/2 rn/ . . , a + 1. 



(vHyV* ^(y) = (- i r( r( a + n + l) ) L " (y)) ^ = 2(n + ^) (2.5) 



where L^(y) denotes the Laguerre polynomial. The factor e E ° T ^ 2 cancels out the formula 
for the distribution function, so the explicit value of E is not needed. 

Knowledge of the Green function allows the p.d.f. 

( (1) (2) (2) (2) \ 

for the event that there are eigenvalues at y^\ . . . , yff when the parameter equals r 1; and 
at y\ , . . . , y N when the parameter is increased to T2 to be calculated. Explicitly 

= ^ I dy[ 0) ■■■j p dyf Po (y{°\ • • • , y®)G{yl°\ • • • , y%\ y?\ ■ ■ ■ , y®; n) 
xG(^ 1) ,---,^ 1) ; 2 /| 2) ,---,^ 2) ;r 2 -r 1 ). (2.6) 

Here po denotes the prescribed initial p.d.f. (in the y variables) and I' denotes the trans- 
formed domain (for convenience this will be omitted from subsequent formulas). Initial 
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p.d.f.'s with orthogonal ((3 = 1) and symplectic (J3 = 4) symmetry are given by the 
functional forms 



., . nf =lv M%)nf<j%--y,i, /? = i 

Po(yi,---,yN) oc ^ v (2.7) 

]lj=i - y j+ N/2) Uj<i \yj -yi\ , P = 4: (N even), 

The parameter dependent distribution function p( n +m) is calculated from the p.d.f. 
according to the formula 

I (i) (1) (2) (2) \ N\ N\ 

P( n +m){yi ,---,vn ,n',vi ,---,y N , T 2> 



(N — n)\ (N — m)\ 

x / dVnli--- / dy [ N / dy^+i--- / dyPp{y{ 1 \---,yP;T 1 ;y[ 2 \---,yP;r 2 ) {21 



In a recent work [|7j we have shown that for initial conditions ( |2.7|) , and G specified 
by ( |2.2|) with general weight function w(y) and corresponding orthogonal polynomials 
{Pk(y)}k=o,i,... , the quantity (|2.8|) can be expressed as a (n + m)-dimensional quaternion 
determinant Tdet (this quantity, which was introduced into random matrix theory by 
Dyson 0, is explicitly defined in e.g. [[|). 

The entries of the quaternion determinant depend on the functions 

S(x,y;r x ,T y ), I(x,y,T x ,T y ), D(x,y;T x ,r y ). 

These functions in turn depend on the quantity F(x, y \ r), which is defined for the (3 = 1 
initial conditions in (|2.7|) by 



F(x, y;r) = j dz j dz(g(x, z; r)g(y, z'\ r) - g(y, z; r)g(x, z'; r)) , (2.9) 



while for the (3 = 4 initial conditions in ( 2.7|) it is specified by 

r d d 

F(x, y;r) = J dz{g(x, z\ r)—g(y, z; r) - g(y, z; r)—g(x, z\ r)}. (2.10) 

They also depend on a family of skew-orthogonal monic polynomials Rk{x; r) of degree 
k. The terminology skew-orthogonal means that they satisfy the skew orthogonality 
relations 

< R 2 m{-;i~), iWi(-5 r) >= - < R 2n +i{-] r), i? 2 m(-; r) >= r m (r)5 min , 

< R2m(-;T),R 2n (-;T) >= 0, < R 2m+1 (-;r),R 2n+1 (-;r) >= (2.11) 



where 

< 



f(x),g(y) >= \fdxj dy^w(x)w(y)F(y, x; r) (f(y)g(x) - f{x)g(y)) . (2.12) 



The polynomials -Rfe(x; r) can be deduced from their static counterparts |10|, Suppose 
at r = we write 

n 

R n (x; 0) = a njCj(x), a nn = 1, (2.13) 

3=0 
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where the Cj(x) are the monic version of the orthogonal polynomials Pj(x) in (|2.3|) . Then 
it is straightforward to verify that 

R n (x; r)e^ T = £ a n3 C 3 (x)e^ T (2.14) 

j=0 

and 

r fe ( r ) = r fc (0)e" (72fe+72fe+l)r (2.15) 
Still another function involved is 



$ fc (x;r) := J F(y, x; r)^f^y)R k {y r) dy. (2.16) 
Now, if we invert fl2.13|) by writing 



C n (x)e^ T = J2Pn 3 e^ T R ] (x;r), (3 nn = 1. (2.17) 

3=0 



then $fe(x; r) can be expanded according to 



/ 00 C (r)e~^ T 



i/=2fc-2 

oo 



2fc-2; 



$ 2fe - 2 (a;;r)e^-^ = ^(^^(O) £ — & zfc-i- (2-18) 

y=2fe-l 

where ft,„ is the normalization 

w(x)Cj(x)Ck(x) dx = hj5j t k- (2.19) 

With this notation, we can now define the functions S, I and D. We have, assuming 
N even, 

S(x,y;r x ,r y ) := yjw(y) 

N/2-1 e ~t2lTx+"t2l+lTy 

x E 7m {®2i(x;T x )R 2 i +1 (y,T y ) -<S> 2 i +1 (x;T x )R 2l (y,T y )}, (2.20) 



1=0 



rj(0) 



I(x,y;T x ,T y ) 

: = E ^ ^(xj^^d/jfj-fwl^^ad/jr)} (2.21) 

l=N/2 Tl ^) 



D(x,y;r x ,T y ) := ^Jw(x)w(y) 

X E {RAx;T)R 2l+1 (yT)-R 2l+1 (x;T)R 2l (yr)} (2.22) 



z=o 



n(o) 



and in terms of S we define 

S(x,y;r x ,r y ) := S(x,y;r x ,r y ) - g(x,y;r x - r y ), r x > r y . 
The formula for P( n+m ) can now be stated. It reads [[7] 



(2.23) 



P(n+m){x\ Xn'', T~i, x\ ' , ■ ■ ■ , Xjf] T 2 ) 

[fx(xf\x^;Ti, Ti)]nxn Lfe (4^ , xj? I r l> r 2)] 

[/f xl 2) ; n, r 2 )]^ xm [A(xf \ xi 2) ; r 2 , r 2 )] 



nxm 
mxm 



(2.24) 



where 



fi(x,y;r,r) :-- 
f2(x,V,T x ,T v ) : = 



S(x,y;r,r) I(x,y;r,r) 
D(x,y;T,r) s(v,x;t,t) 

S(x,y,T x ,T y ) I(x,y;r X: T y ) 
D(x,y,T x ,T y ) S(y,x;T y ,T x ) 



and 



S(y,x;r y ,T x ) -I(x,y,T x ,r y ) 
-D(x,y,T x ,r y ) S(x,y;r x ,T y ) 

If we denote the matrix in ( p. 24 ) by X, and note from the definitions ( 2.21| ) and (|2.22 ) 
that 

I(x, y; t x , T y ) = -I(y, x; r y , r x ), D(x, y; t x , t v ) = -D(y, x; T y , r x ) (2.25) 
then we see that ZX where 



Z = 1 



(ra+m) 



-1 

1 



is an even dimensional antisymmetric matrix. Its Pfaffian is thus well defined, and in fact 



TdetX = Pf(ZX). 
This formula can be used to compute the formula 



(2.26) 



2.1 Alternative formulas and inter-relationships 



As first noted in [|J, substitution of (|2.18|) in ( p.20|) and use of ( P.17|) implies the decom- 
position 



S(x,y;r x ,Ty) = Si(x,y,T x ,T v ) + S 2 {x,y,T x ,T y ) 



where 



.N-l 



Sx{x,y,T x ,T y ) = ^Jw(x)w(y) ^ 



C v (x)C v (y)z~ lv{Tx - Ty) 



(2.27) 



(2.28) 



oo N-l Q( x \ e -lur x 

S 2 (x,y,T x ,T y ) = s/w(x)w(y)J2 



u=N k=0 



h, 



■PvkRk{v,Ty)e™ (2.29) 
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Since < 71 < 72 < • ■ • (recall (|2.4|) and (|2.5|) ) we see that for iV — ► 00 and with all 
other parameters fixed, S2 — > while 

Si (x, y; r x , T y ) -> g (x, y; r x - r v ) , 

assuming t x — r y > 0. The use of this result is that it shows 

I(x,y,T x ,T y ) = J F(y',y;Ty)(g(x : y';r x -Ty) - S(x,y';T x ,T y )^j dy (2.30) 



valid for t x — r y > (for t x — r y < we can first use the first equation in ( |2.25|) , then 
apply O). 

We would also like to relate D to S. Apparently this can be done naturally only when 
t x = 0. Using the fact that g(x, y\ 0) = 5(x — y) we see from (|2.9| ) and ( p.!0|) that 



Recalling (|2.16|) , these formulas show 



d 

dx 



$,(x;0) = 2^{x)R k {x;0) 



$K^0) = -2^{^(x)R k (x;0)) (2.32) 



dx 

in the two cases respectively. Thus for (3 = 1 initial conditions 

d 

dx 

while for (3 = 4 initial conditions 



S{x,y;0,T v ) = -2D(x,y;0,T v ) (2.33) 



f 5(x', y; 0, r v ) da/ = 2D{x, y; 0, r v ). (2.34) 
Jy 

3 The Laguerre ensemble 

For the Laguerre ensemble the weight function w(x) and constants 7^ are given by ( p.5| ), 
while the monic orthogonal polynomials and corresponding normalization have the explicit 
form 

C n (x) = n\(-l) n L a n (x), h n = T(n + l)T(n + a + l). (3.1) 

This applies independent of the particular initial condition (orthogonal symmetry {(3=1) 
or symplectic symmetry {(3 = 4)). However, the skew orthogonal polynomials Rk{x;r) 
and the quantities j3jk depend on the initial condition. 
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3.1 (3 = 1 initial conditions 

In this case the monic skew orthogonal polynomials for r = are known to be 



d lm (-1) 

R 2m {x) = -(2m)\-L a 2m+1 (x) = (2m)\J2^-C 3 
ax j_ j. 



R2m+i(x) = -(2m+l)\L a 2m+1 (x)-(2m)\(a + 2m + l)^L% m (x) 

= C 2m+1 (x) + {2m)\(a + 2m+l) £ ^-C^x). (3.2) 

with corresponding normalization 

r n (0) = 4T(2n + l)T(a + 2n + 2). (3.3) 

These equations can readily be inverted. We find 

„ ,s „ / n /o m m! v^ 1 (j + a/2 + l)!(2j + 2)! 
C 2m (x) = iW*) + (2m + a)\ 1 -^ m E (J + ; )!(2J } + \ + 2) / 

X ((27Tiy^ +l( ^"(2^!^ (x; 

C 2 m +1 (x) = iWi(s) + (2m + a + 1)! ^ £ U + ^ + ^ 

(m + a/2)!^ (j + l)!(2j + a + 2)! 

X ((27T1)!^ +1(X) -(2^^ (X) ) (3 ' 4) 

(direct sustitution of (|3.2|) into the right hand sides verifies the validity of these formulas). 
According to the definition ( |2.17| ) with r = 0, we read off from ( |3.4| ) that 

$2m 2m = /^2m+l 2m+l = 1; P2m+12m = 

3 - u-iav 1 [( " /2)]! ^ + a / 2 + 1 ) ! ^ + 2 ) ! rssi 

Pu2i ~ {U + a) -(2j)l([^)+a/2)\ + (2j + a + 2)\ ^ 

where v > 2j + 1. A most significant feature of these formulas is that the dependence on 
p and / in f3 p i separates. In particular, for n > N — 1 and N even we have 

(3 P i = P p n-iPn-u, 0<1<N-1, l^N-2. (3.7) 



Substituting this in the formula (|2.29| ) for S 2 , and making use of the final formula in (|3.5| 
in the case (v, 2j + 1) = (p, N — 1) gives 

C p (x) 



S 2 (x,y;T x ,T y ) = <Jw(x)w(y) ~t — e lpTx P P N-\ 

p=N n P 



N-l 

X ' 

1=0 



1 v — ± 

( £ p N -uHv\ Ty)e llTv ~(N- l)R N -2(y, r y )e^-^). (3.8) 
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Now the sum over I can be computed using (|2.17 ), while the sum over p is given by ( p. 18 
Thus we have the explicit evaluation 



S 2 (x,y,T x ,Ty) = yjw(x){ — 77v\ \/w(X; 

x (C N ^{y)e^-^ - (N - l)i2*_ 2 (y; r,)e^). (3.9) 

In Appendix A we will show that in the limit t x , r y — > this expression reduces down to 
the form of 5*2 recently derived by Widom || . 

3.1.1 The thermodynamic limit 

In the neighbourhood of the spectrum edge at the origin, we will show that after intro- 
ducing the scaled coordinates and parameters X and t according to 

x = — r=— (3.10) 

(the numerical factors 1/4 and 1/2 are chosen for convenience), the distribution functions 
have a well defined large N limit. In particular, the scaled distribution function ( |2.8| ), to 
be denoted pfe? m ), is given by 

hard (x (l) vil). f {l). X m X^-t (2) ) 

= lim { ±y^ P{n+ j^...^^ HA (3,!) 

n^ooK^N) IJ{ - n+m A^x' '4A'4A' 4A ' ' 4A ' V ; 



From the formula ( 2.26|) for Tdet it follows that p\n+m) ls gi ven by the formula Q2.24 ) with 



the functions S, S, I, D replaced by their scaled counterparts 



S^ d (X,Y;t x ,t Y ) := lim — s(— ,— ;—,—) (3.12) 

S^(X,Y;t x ,t Y ) := (3.13) 

I^(X,Y;t x ,t Y ) := & J(^, (3 - 14) 

D h ™ d (X,Y;t x ,t Y ) := lim f— )V— , — ; ^, (3.15) 

where the subscripts 1 indicate (3 = 1 initial conditions. 

Note that it was permissable to multiply -^1 by 4A and divide j^D by 4A in the 

above definitions because of an invariance property of the Tdet formula (|2.24j) . Thus in 
general each 2x2 block Qjj. can be replaced by the 2x2 block 

A^Q jk A. (3.16) 

Choosing 



V4A 

then transforms j^I and -^D as required. 
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Asymptotics of S(x, y; r x , r y ) 

Let us first compute the asymptotic behaviour of $2^-2 r )- This is required in the 
computation of S2 (with 2k = N) , and also I. From the expansion ( 2.18|) and the explicit 
formulas (|3.1|) and fl3.5|) for the quantities C v (x), h v and /3 U 2k-2 therein we have 



L 



2k-l 



-72fc-lT 



<™W ! ^ " -i; +lW e— )).(3.18) 



+ 



E 



(2A;-l)!(2fc + a)!fc!^(i/ + a/2)! 

We now proceed to estimate the summand for large v. 
From Stirling's formula 

v\ 1 



(i/ + o/2)! 1W 2 ' 
For the second factor in the summand we write 

L a 2[/ (x)e-^ T - L a 2u+1 (x)e-^ T 

= -e-^ T (L a 2u+1 (x) - L a 2u (x) - L a 2u (x)(e- T ^-^ - 1 

Using the facts that 



gT(7 2 „-72„+i) _ g-2r 



L 2u+ i(x) — L%,,(x) — L 



2iA 



J 2u+l( X ) 5 



and the asymptotic formula 

^)L a n {x) n ~ oo n a ' 2 J a (2^) 
where J a (z) denotes the Bessel function, we see that 



(3.19) 



(3.20) 



(3.21) 



w\x) 



{y + a/2)\ 



L a 2v {x)e-^ T ~ ja 



L 2vJr i{x)e 



x=X/4N 
r=t/2N 



~ -v~ a l\- 



u:=2v/N. (3.22) 



(The equality is obtained by making use of a suitable differentiation formula for the Bessel 
function). 

Substituting (|3.22 ) in ( 3.18Q gives the leading order behaviour of the sum over v as a 
Riemann integral. Furthermore, use of ( |3.21[ ) gives 



'w(x) 



^2k-l( X ) 

(2k-l + a)\ 



"72fc-lT 



jy-o/2 



x=X/iN (2k 

T=t/2N 



— s- a l 2 J a {yiX) e 



St 



(3.23) 
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s := 2k/N, so in total we have 



2k-2 



\4N 2NJ 



,t2k- 2 t/2N 



r - ( °)(S^( 1 - 2sW2 ) J «(^)^ 

g— ut ^ 
-J a (VuX) du 



(3.24) 



The expansions ( |3.23| ) and ( |3.24j ) with 2k = N and thus s = 1 suffice to specify the 
^-dependent factor in (|3.9|). Thus we have 

■C N ^{x)e-^ N -^ 



jy-o/2 



AT-l 



x=X/4JV 
T iE =t x /2JV 



2(JV- 



£i / poo o — LLL X \ 



(3.25) 



Next we will compute the asymptotic form of R 2 j{y; T y ). In the case 2j = N — 2 this 
is required for the ^/-dependent factor in (|3.9|) . From (|2.14|) , ( |3.2|) and ( |3.1|) we see that 



2j 



i2 ! y(y;r)e'W = (2j)!E^(f) ^ ■ 



(3.26) 



i=0 



The asymptotic form of this expression is deduced by using ( |3.21|) to approximate the 
summand, which leads to the result 



y=Y/4N 
Ty=ty /2JV 



(2j)!iV 1+a/2 f v a/2 J a (VvY) 



dv . 



(3.27) 



where here s = 2j/N. Setting s = 1 and making further use of ( |3.23|) we see that 
V^(y)(C W -i(y)e^- lT » - (N-l)R N _ 2 (y;r y )e^-^) 

~ -(iV - l)!iV 1+a/2 jf* v a/2 J a (v^F) cfo . (3.28) 
Combining fl3.25| ) and ( |3.28j ) in the formula ( |3.9| ) for S 2 gives 



x U\ a ' 2 J a (VvY) e vtY dv^j . 



(3.29) 



The asymptotics of Si can be determined by using ( |3.21| ) and Stirling's formula to 
approximate the integrand. This gives 



\4N' 4JV' 2N' 2NJ 



N jf* J a (VuX) J a (VvY) e - u{tx - tY) du . 



(3.30) 



Since S is defined as the sum of Si and S 2 , its asymptotic form is given by the sum of 
(|3~29p and QOOl) . Substituting into fl3T2D gives 

S^ d (X,Y;t x ,t Y ) = - t J a (V^X)J a (VuY)e- u ^-^ du 

4 jo 

1 , , a roc p—utx . i-l , 

+ -(-e- t *J a {^fX) + - J a (VuTx)du) dvJ a (VvY)e tYV v a / 2 .(3.31) 

8 v 2 Ji u ' Jo 
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Asymptotics of I{x,y;r x ,T y ) 

So as to induce a cancellation at leading order, we first perform some minor manipulation 
to the definition (|2.21| ), by way of subtracting and adding (2l+i)e y2l ( T * +T ^$ 2 i(x', i~ x )$ 2 i(y, 
to the summand, so that it reads 

oo 



r„ 



1 



I{x, y] T x ,r y ) = Y, -j^{e l2lTx $2i(x-,T x )(e^ T y$ 2l+1 (y-,T y 

l=N/2 r 'V U J 

-(2l + l)e^ T ^ 2l (y;r y ))-(x^y)} . 



(3.32) 



The asymptotics of e l2lTx $2i(%] t x ) can be read off from ( |3.24j ). For the ^-dependent factor 
of the first term of the summand we note from the definitions Q2.18 ) and the explicit 
formulas (B.ll) and (B.5I) that 



e 72i+lT $ 2 /+i(y; t) - {21 + l)e^ T $ 2l (y; r) 

= V«(y)n(o) [(2i + 1) {2l + 1 + a)] 



Thus, using ( |3.3| ) and the asymptotic formula ( |3.21| ) we have 
(e™+^$ 2m (y; r) - {21 + l)e»' T $ a fe; r) 

~ 4(2/ + l)!AT- 1+a / 2 | A (V s * J a (v 7 ^)) - as- 1+a / 2 J a e~ s *| , s = 2Z/JV. 



y=Y/AN 

T=i/2JV 



Substituting this result, together with ( |3.24j) (after replacing k by Z + 1) in (|3.32|) shows 
that 



I^(X, Y; t x , t Y ) =fds\(l- J a (v^x) e"** - | 



oo g—utx 



du 



1^ (V^J. (VsP)) - as" 1+a/2 J a (ViP) e~ s ' y j - (X <-> F$3.33) 



Asymptotics of g(x,y; r 

According to ( |2.3| ) and 



^(x, y; r) = (xy) a/ V ( ^ )/2 V — ^ . , . 



L J a (x)L J a ( 2 /)e- 2 ^' + ( a+1 )/ 2 ) r 



(3.34) 



This summation can be evaluated in closed form 11 



giving 



5 



(ar, y; T ) = e -(*+y)/2 e -(2o+i)r (i _ e -2r) 



-1 



-2t 



exp 



l_ e -2r^ ^ «y l _ e -2r 

where J a denotes the Bessel function of pure imaginary argument. Consequently 
(X Y tx-ty 



y 



V4iV'4iV' 2N 



VXY 

tx—ty V tx—tyJ \2{tx—tY) 



exp ( — — ) /„ 



(3.35) 
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Alternatively, proceeding as in the derivation of ( |3.30|) we can deduce from (|3.34j) that 

X Y t X -ty. s N 



AN' AN' 2N 



r J a {\f^X)J a {VuY)e~ u{tx - tY) du. 
Jo 



(3.36) 



The immediate merit of this representation is that when substituted in ( |2.23j ), partial 
cancellation occurs with ( |3.30 ). Adding the result to ( |3.29| ) and substituting in ( 3.13 ) 
gives 



S^ d (X,Y;t x ,t Y 
1 



1 

4 Ji 



J a (VuX)J a {VvY)e- u{tx - tY) du 



e- tx J a (Vx) + - J a (VuX)du)+ dvJ a (VvY)e tYV v a/2 (3.37) 

2 Jx u ' Jo 



Asymptotics of D(x,y,T x ,T y ) 

According to (|2.22j ) D is defined in terms of {Rj(x; r)}. Now R 2 j(y;r) is specified by 
( |3.26j ). Deriving the analogous formula for R 2 j + i(x;r), which follows from (|3.2|) , (|2.14|) 
and (|3.1|) , we see that 



R 



2 j+ i{x;T)e 



72j + l1 



-(2j + l)\L a 2j+1 (x)e^ T - (2j)\(2j + a + l)L%(x)e^ 
+{2j + a + l)R 2j {x; r)e^ T . 



Hence 



(e^ T *R 2j+1 (x; T x )e™R 2j {y; r y ) - {x <-> y)) = 

(-(2j + l)\L a 2j+1 (x)e^+^ - (2j)\(2j + a + l)L a 2j (x)e^) R 2j (y; T y )e™ - (x <- y) 

Now the asymptotic form of R 2 j(y; r y ) is specified by ( p.27| ) while we see from (|3.21|) that 



(-(2j + l)\L% +l {x)e^ T * - (2j)\(2j + a + l)L%{x)e^ 
~ -2(2j + l)!e st *(As) a/2 J a (J^x) , s = 2j/N. 



x=X/4N 
T X =t x /2N 



Making use also of the explicit form ( |3.3|) of r n (0) and using the definition ( |3.15| ) we thus 
find 

D^ d (X,Y;t x ,t Y ) = -A~ 3 

x y dse stx s- a/2 J a (V^X) J° dve vtY v a/2 J a (VvY) - {X <-> F)| . (3.38) 

3.2 (3 = 4 initial conditions 

The method of evaluation of the quantities S, I, and D for the Laguerre ensemble in the 
case of initial conditions with symplectic symmetry (/3 = 4) closely parallels that just 
given for (3 = 1 initial conditions. First, we note from previous work || that the monic 
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skew orthogonal polynomials at r = have the explicit form 

«*<*> = ^•■fa" + (- 1 ^'S^ ( t+ (.li)/2 )l ^ 1 w 

= 2*JI + (« - D/2)! £ ; , 22l(; + ( 1 a _ 1)/2) , (0,(«) + 2K^ l( ,)) , 

*aj+iM = -(2j + l)!L-i 1 (x) 

= C 2j -n(x) + (2j + l)C 2j (x), (3.39) 
with corresponding normalization 

r n (0) = T(2n + 2)T(a + 2n + l) . (3.40) 

The inversion of these equations gives that the quantities (3 n j in ( |2.17| ) have the explicit 
values 

/?2m2m = /?2m+l 2m+l = 1 , ftm+12m = —(2m + 1) , 

(3 v2 j+\ = (—l) u+1 ^. — ^ , /?j,2j = aj3 u2 j+i , (3-41) 

for v > 2j + 1. Thus the factorization formula (|3.7| ) again holds, which when substituted 
in (|2.29|) and after making use of the explicit value of (3n-in-2 and (3 p n-i, and ( |2.17|) , 
implies 

/ 00 C (t) 

S 2 (x,y;r x ,T y ) = ^Jw(x)w(y) £ -f-L e">- (3 pN ^ {C N - X {y)e™-^ 

p=N n P 

+ (a + N-l)R N ^(y;T y )e^ T y} . 

The sum over p is, after minor manipulation, formally identical to the sum over p in ( |3.8| ) . 
so it can be summed using fl2.18|) . Thus S 2 has the explicit form 



\ ?"jV/2-l(U) / 

x (CW-i(y)e w - lT » + (a + iV - l)ifo_ 2 (y; rje^ 21 ') . (3.42) 

As with (|3.9|) , in the limit r x , r y — > this reduces down to the form of S 2 derived by Widom 
|| in the static /? = 4 theory. The details of the verification are given in Appendix A. 



3.2.1 The thermodynamic limit 



Although the scaled distribution function is still given by ( |3.11| ), the required scaling 
of the elements of the quaternion determinant ( |2.24j ) differs slightly from the formulas 
( 3.12 )— ( 3.15 ) appropriate with (3 = 1 initial conditions. In particular Q3.14 ) and (|3.15 ) 
should now read 

lim f — 1 l(——;^,JL) (3.43) 



I^ d (X,Y;t x ,t Y ) 



A'- 



D^ d (X,Y;t x ,t Y ) 



(3.44) 



' V4iV' AN 2N 2NJ 
AN' 2N' 2N) 

while (|3.12|) and (|3.13|) remain the same, except that the subscripts 1 on the right hand 
sides should be replaced by 4. 



lim D(—. 
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Asymptotics of S(x, y; r x , r y ) 

From the expansion ( |2.18| ) and the explicit formulas ( |3.1| ) and ( |3.41| ) we have 



r k -i(0)y/w(x) 
(2k- 1)1 



oo j a 



v=2k-l 



Applying the asymptotic expansion (|3.21|) and Stirling's formula shows that 



$2fc-2 



UiV' 2NJ 



72fc-2*/2iV _ T-fc-l(O) jyl-o/2 



(2Jfe-l) 



y°° rf rr -«/2 e -^ j a ( v ^x) ) (3.45) 



where w := 2k/N. Setting it = 1 gives the asymptotic form of $N-2{%)T x )/rN/2-i(0) as 
required in (|3.42|) , and comparison with (|3.23| ) shows that this term dominates the one 
involving Cn-i(x). Thus here the formula ( |3.25|) applies with the right hand side replaced 
by 



Arl-a/2 „oo 

/ drr~ a/2 e~ rt T ( vrX 

For the y-dependent factor in ( |3.42|) , we note from (|2.14|) , ( |3.1| ), and ( |3.26| ) that 

2^! 0- + (a - l)/2). ± ( i5lfe)e w - 

Use of (|3.21|) and Stirling's formula gives for the corresponding asymptotic form 



(3.46) 



y=Y/4N 
T=t/2N 



I" \~a/2d Lut^j^ Lf^\\ du ? g . = 2 j/jV. (3.47) 
1/ (XIX 



W2-1 

-l+a/2 



Adding this with s = 1 to the asymptotic form of \J w(y)C n _i(y)e lN ~ lTy / (N — 1)! deduced 
from ( |3.23| ), and multiplying the result by ( |3.46| ) gives that 



2 UiV'4iV'2iV'2iVy Ji ^{yr^j 

x (-J a (VF) e* + i jf 1 M -/ 2 ^ (e«^/ 2 J a (v^l 7 )) du) . (3.48) 

The asymptotic form of Si can be read off from ( |3.3(J| ) since according to its definition 
( |2.28j ), its value is the same for (3 = 1 and 4. The asymptotic form of S itself is thus given 
by the sum of (|3.30|) and fl3.48| ). Thus we have 

S^ d (X,Y;t x ,t Y ) = 

l r 1 



A Jo 
xf- 



J a (VuX)J a (VvY)e- u{tx - tY) du-- drr- a/2 e- rtx J a (V^X) 



4 Ji 



J a (VY)e tY + - f\- a/2 ^(e utY u a/2 J a (VuY))du 
2 Jo du K ' 



(3.49) 
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Asymptotics of S(x, y; r x , r y ) 



The single particle Green function g as specified by ( |2.3| ) is independent of the initial 
condition, so me can use the formula (|3.36| ). Dividing by AN and subtracting from (|3.49 ) 
shows 

S^ d (X,Y;t x ,t Y ) 

1 roo 1 roc 

= --/ J a {V^X)J a {^vY)e- u{tx - tY) du- - / drr~ a/2 e~ rtx J a (V^X) 
4 Ji 4 Ji 

- J a (VY)e tY + 1 j\- a l 2 ^-{e ut -u a l 2 J a {V^))dv)j (3.50) 



Asymptotics of I(x, y; t x , T y 



The definition ( |2.21[ ) gives that I is defined in terms of {§i(x; r)}. The asymptotic form of 
&2k-2(%'i T ) is given by ( [3.45 ). Instead of deriving an analogous expression for $ 2 fc-i(z/; T )i 
we note from (gTg ), ( pOT] ), and the formula ( ggg ) for $ 2fc _ 2 that 



^afc-i(y;T)e 



72fc-lT 



Jw(y)r k -i(0) , , 

( 2A . _ 2 + a) ; (^V 2 (yK^- 2T + ^-iCyJc-**-^) - a^Cy; r) . (3.51) 



Thus 



e™ T *$ 2 ,(x; T^eW+^a-fiCs/; r„) - (re <-► y) 

/w(y)r/(0) 



(2/ + a)! 



j (^(y)e-^^ + LVi(y)e- 72!+lT «) - (s «-> y). 



(3.52) 



According to ( |3.21| ) and ( |3.40| ), the product of the second and third factors in the first 
term have the same leading asymptotic behaviour 



- (21 + 1)! A(vN) a/2 J a (VvY) e~ VT , v = 21/N . 



(3.53) 



Multiplying this by ( |3.45| ) (with k \— > I + 1) and substituting in (|2.21 ) gives that the 
desired asymptotic form is 



/ 4 hard (X,F;t x ,t y ) 



I {/°° dT r ~ a/2e ~ Hx J « ( v ^) } ^ J " ( v ^ 7 ) e ~" ty - (X <-> F).(3.54) 



Asymptotics of D(x, y\ t x , r y ) 



Here we require the asymptotic form of {Ri(x; r)}. For I even this is given by (|3.47 ). On 
the other hand, for I odd we see from ([2 . 1 4j) , ( |3.1| ) and ( |3.39| ) that 



R2 j+ i(y\ r)e^ = -(2j + 1)! (L a 2j+1 (y)e^ - L%(y)e™ 



(3.55) 
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Use of (ET23D shows 

fMjjj Rq+ifa r)e^+- y=Y/AN ~ -(2j + l)!^^ 1 ^ (eW a , (3.56) 

=t/2N 



s := 2j/N. Multiplying this by (2j + 1) ! times ( 3.47] ) (with Y" replaced by X) and dividing 
by ( p.40p we thus have 

D 4 hard (X, Y; t x , ty) \ f Q ds s- a ' 2 1 jf* u- a ' 2 ^ (e ut *u a ' 2 J a (v^x)) du\ 

x Ts i eStYsa/2j - (^)) - ( x ~ r ) • ( 3 - 57 ) 
4 Gaussian ensemble 



For the Gaussian ensemble, we see from ( 2.4Q that the corresponding monic orthogonal 
polynomials and normalization are 

C„(x) = 2~"i/ n (y), /i n = 2-" 7 r 1 / 4 n!. (4.1) 

The skew orthogonal polynomials for both (3 = 1 and (3 = 4 initial conditions at r = are 
known in terms of the C n (x) (see e.g. [P^j). Using these expansions, the general formula 



( 2.14 ) then gives 



R 2m {x;r)e^ T = C 2m [x)e^ T (4.2) 

R 2m+ i(x;r)e^ T = C 2m+1 {x)e^ T - mC 2m ^(x)e^ T (4.3) 

r m (0) = 2- 2m+1 v^Fr(2m+l) (4.4) 

for (3 = 1 initial conditions, and 

R 2m {x;r)e^ = m!£e^^p 

^ 2m+1 (x;r)e^ +1T = C 2m+ i(*)e 72m+1T (4.5) 

r m (0) = 2- 2m ^T(2m + 2) (4.6) 

for /5 = 4 initial conditions. Inverting these formulas for {Rj(x)} in terms of {Ci(x)} 
according to (|2.17|) gives 



f m\/l\, j = 21 + 1, 1 = 0,1,... ,m , s 

fhmj-dsmj, /Wii-| Qj j = 2Z, Z = 0,l,...,m (4 ' 7j 

for (3 = 1, while for /3 = 4 we find 

1, j = 2m 

/32m j = "( -»™, j = 2m - 2 /?2m+lj = <WflJ. (4.8) 

0, otherwise 

Here we will use these explicit values in the general formulas of Section 2 to compute 
the matrix elements in ( [2.24Q . In particular, we will give their asymptotic form in the 
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neighbourhood of the soft edge, which is specified by introducing the scaled variables X 
and t according to 



x = V2N + 



X 



t 



(4.9) 



21/2 jyi/6 ' jyi/3 ' 

The corresponding scaled distribution function ( |2.8|) , to be denoted Pu^^, is given by 

soft (x (l) x (l).(l). x (2) ¥-(2).f(2h 



/ 1 \ n+m / I 

= ^SW^) P(n+m)(V2N + 



(1) 



21/2JVV6' 



, tW , X) 



(2) 



x 2iV + 



X^ t (2) 



(4.10) 



21/2 N 1 / 6 ' N l/3 ' v *" ' 1 2 V2j\rl/6'"'' 2 1 / 2 A^ 1 /6' ]\Tl/3, 

Analogous to (|3.12|) - (|3.15|) , the matrix elements of p^f+ m ) are computed from the formulas 

Sl° (X, Y;tx,ty) 



e -N 2 ^(t Y -t x ) X 

lim 77T 777 S[v2N H 777 7777. 



y 



S?*(X,Y;t x ,t Y ) 



e -N 2 / 3 (t Y -t x ) _ 



2l/2jyl/6 ; 



lim e 

iV— >oo 



I(V2N + — 
\ 01 



e -iV 2 / 3 (* y +tx) 

lim -77 

X^oo 2A^ 1 /3 



21/2 jyl/B' 



2V2AT1/6 


' jyi/3' jyi/3 


y 


ty 


21/2JVV6 


' jyi/s ' jyi/3 


y 




2l/2jyl/6' 


jyi/3' JVI/3J 


y 


tx ty 



(4.11) 



) (4-12) 



(4.13) 



2 i/2jyi/6 



21/2JV1/6' j\Tl/3' jyl/3 



Here the scale factors in addition to l/2 1 ^ 2 N 1 ^ 6 do not effect the value of ( |2.24 ) because 
they result from a transformation of the form ( |3.16| ). 



4.1 (3 = 1 initial conditions 

We see from (|4.7| ) that for N even and p > N 

(3 p i = PpN-i Pn-u ■ (4.15) 
Substituting in the formula (2.29) for S2 and using ( J2.17|) and ( [2.18D we thus have 



S 2 (x,y,T x ,T y ) = ^w(y)( 



rjv/ 2 -i(0) 
xC N ^(y)e^-^. 



C M -i(i)e" w - lTl 



h 



N-l 



(4.16) 
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To compute the scaled limit we make use of the asymptotic expansion [13| 

e~ x2/2 H n (x) = vr- 3 / 4 2"/ 2+1 / 4 (n!) 1 / 2 n- 1 / 12 (vrAi(- M ) + 0(n- 2/3 )) (4.17) 
where x = (277,) 1 / 2 — u/2 l / 2 n 1 ^ . This gives 



f^ C »M yHm v> + r,. W ~ * 1/4 2-^y^(N - l)^N-^M(Y)- (4.18) 
Also, from ( gig) , (pj) and 

r-jv/2-i(0) v /ijv-i 

" (iV/2-1)! 6 p ^ o 0F(iV + 2p + l)!^ +2p+l(x)e 
Noting that for large N 

(N/2+p)\ ((iV-1)!) 1 / 2 (p+1) 
((JV + 2p+l)!)V2 (JV/2-1)! 

and making use of (|4.18 ) with n — N + 2p + 1, — it ~ X - 2p/N 1/3 shows 

X t x 



Hence 



((N-iy.)^f((2N)^ + 21/2N1/e , N1/3 

^ e -«W, 2 (JV-l)/2 2 l/4 ;r -l/4 JV -l/l 2 ^! /" Ai(Jf _ (419) 

2 jo 

+ W^W' < 2iV ) 1/2 + 2^ i^' If*) 
~ e 7v2/3( * y -* x) 2- 1/2 X 1/6 Ai(r) f° Ai(A - u)e- wtx tZv. (4.20) 



Now consider Si as defined by ( 2.28|) . Substituting (|4~T|) we see that the asymptotics 
can be determined by writing the sum so that the sum index p is replaced by N — p, and 
then using Stirling's formula and ( 4.18 ). We thus find 

5 1 ((2iV) 1 / 2 + 2l/2]Vl/6 , (2A0 1 / 2 + 2l/2jVl/6 ; — j^, -^3) 



_ e 7v 2 /3 (iy -^) 2 i/ 2iV i/6 / Ai / X + ^) Ai (y + ^W^"**) du. (4.21) 

Jo 

Adding together ( |4.20|) and (|4.21| ) and substituting for 5 in Q4.11Q we thus have 

poo I /-co 

^(X.yjt^ty) = / Ai(X + v)Ai(Y + v)e- v( - tY ~ tx) dv + -Ai(Y) / Ai(X— ^e^* cfo. 
jo 2 Jo 

(4.22) 
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Asymptotics of I(x,y;T x ,r y ) 

We know from ( |2.21| ) that / is defined in terms of {$;(■; r)}. The asymptotics of <& 2 k~2 is 
calculated in an analogous fashion to the expansion ( |4.19| ). Thus after noting 

= r N/2+ k-i(0) Nr y> H 2m+N+1 (x)(N/2 + m)! (2m+3/2)r 
(iV/2 + A;-l)! j£_ x 7r 1 /2(iv + 2m + l)! 

we find 

<$> N+2k M2N) 1 ' 2 + X ■ J^) e ^ +2k - 2 t x /N^ „ r N/2+ k-i(0) (N/2-l)\ 
WN+2k-2\K"v ) t- 2 i/ 2JV -i/6' ivi/3; e ((7V-l)!)V2(jv/2 + fc-l)! 

Arl/3 /»oo 

xe -^/ 3 *X2^-l)/2 2 V4 7r -l/4 iV -l/12^_ ^ ^ Ai(x _ w)e -, ix dv ^, 23 ) 

Consider now $ 2 fc-i- We note from (|2.18|) , Q and ([O]) that 



HN+2k-2{y) 



so use of (|4.17| ) gives 



xe -2 te /J/l/S ;r _ 1/ 4 2 ( J , +tt _ 2)/2+ l / 4 ((JV + 2j( . _ 2)1)1/2^-1/12 Ai( y _ (4 ^ 4) 



Multiplying (|4.23|) and (|4.24j) together, substituting in ( [2.21| ) and making use of the 



explicit value of r m (0) from (|4.4j) , we then find 

If° h (X,Y;tx,tY) = ds Ai{Y - s)e~ stY Ai(X - v)e~ vt "dv - (X Y) (4.25) 

JO Js 

Asymptotics of D(x, y\ t x , r y ) 



Here the asymptotics of {R m {-\ r)} are required. Now, from (|2.14 ) and (|4.4| ) we have 



RN-2i-2(x;r)e^-^ T = 2"( Af - 2 ^ 2 ) H N _ 2l Jx)e^-^ T . 



Use of (|4.171) then gives 

7N-21-2T _ c N 2 / 3 t x 



w{x)R N - 2 i- 2 (x;T)e 



x .=(2Af) 1 /2 +X/2 l/2 jV l/6 ~ e 



r=t x /V 1 /3 
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Furthermore 



Xe -2 ltx /^/a 7r l/4 2 -(JV-2I-a)/2 2 l/4 ( ( JV _ 2l _ 2)1)1/2^-1/12^^ + ^g) 



R 2m+1 (x;r)e^ T = 2^ 2m+1 ^ H 2m+1 (x)e^ T - m2^ 2m ^ H 2m ^(x)e 



72m-lT 



_ e 72m+ir e x 2 /2_^_/ e -x 2 /22-2m^ ( 
rl/r V 



2m(X, 



-m2~ {2rn - l) H 2m ^(x)e r2m - 1T {I - e 2r ), (4.27) 
21 



and use of ( |4.17|) gives for the asymptotics 



^/Hx)R N -2i-i(y-,r)e^- 2 ^ T 



N 2 /H Y 

j/=(2Ar) 1 /2 + y/ 2 i/2Ari/6 ~ e 
T =t Y /N 1 / 3 



X7r l/4 2 -(Ar-2i-2)/2 2 3/4 ((iV _ % _ 2) ;) 1/2^1/12^. (g-^ Ai (y + s) \ ( 4 . 28 ) 

s := 21/N. Multiplying together ( |4.26|) and ( |4.28| ), dividing by r^ N _ 2 i~2)/2 an d summing 
over I (which is of the form of a Riemann integral) we thus obtain 

DT h {X,Y;t x ,t Y ) = -J dse- stx ki{X + s) — (e- stY ki{Y + s))-{X^Y) (4.29) 

Asymptotics of g(x, y; r y — t x ) 
According to (|2.3|) and 



This summation can be evaluated in closed form 11 



giving 



g — t/2 g — 2r 2e _ T 

g(x, y; r) = e -(- 2 +y 2 )/^-—— __ ex p ( - 1 _ e _ 2r (x ~ V?) exp ( ^——^ xy) 
Hence 

q((2N) 1 / 2 + X (2N) 1 ' 2 + Y ■ —) ~ ivV6 -(x-y)V4M(i+y)/2 



Alternatively, proceeding as in the derivation of (|4.21|) , we can deduce from (|4.30| ) that 
9 « 2/V)1/2 + ^- (2Ar)1/2+ V2W^) 



e 



-N>/H 2 i/2 N i/6 r Ai ( X + v ) Ai (y + v yt dv ( 4 31 ) 

■/ — oo 



Dividing by e Ar2/3 *2 1 / 2 A r1 / 6 , replacing t by tx — £y and subtracting from ( 4.22|) shows 



Sf\X,Y-t x ,tY) = ~ [° Ai(X + v)Ai(Y + v)e^ v{tY ~ tx) dv 

J —oo 

1 r°° 

+-Ai(Y) / Ai(X - u)e- wtx du. (4.32) 

2 jo 

4.2 (3 = 4: initial conditions 

Here we want to again compute the scaled distribution (|4.10|) . As for (3 = 4 initial 
conditions in the Laguerre ensemble, we must modify the scale factor in ( 4.13Q and Q4.14 ) 
so that 

J 4 soft (A,F;tx,ty) 

7V 2 / 3 (t x +iy) _ X Y t 

= lim — - — l((2N) 1/2 + -= -,(2N) 1/2 + -= -: ——) (4.33) 
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HF{X,Y-t x ,t Y ) 

= lim e- N2/3 ^ +t ^D((2Ny/ 2 + * - , {2N) l l 2 + -=?—- ; -J^) (4.34) 

Substituting Q4.8| ) in (|2.29|) shows that in terms of {Ci(x)}, S 2 only consists of a single 
term, which reads explicitly 

S 2 (x,y;r x ,T y ) = ^{x^y f^^' ( - y)ifo- a (y; r y )e™-* r : (4.35) 

The explicit form of Rn-2{v', T y )e 1N ~ 2Tv is given by (|4.6| ). This expression however does 
not give a Riemann sum after substituting the asymptotic form of the integrand. Instead, 
we make use of the identity 

J2 ^2- 2l H 2l (x)c~ l = -a(Ac)- m e ax2 f e~ as ' H 2m+1 (s) ds, (4.36) 
i=o 11 J - QC 

c = (1 — a) /a. This can be verified by multiplying both sides by e~ ax2 and differentiating 
(we also verify that both sides have the same x — > oo behaviour). Choosing c = e~ 2r and 
comparing with (|4.6j ) we see that 

R 2m (y;r)e^ T = --^^(2e~T 2m e y2/il+e ' 2T) r e- s2 ^ 1+ ^H 2m+1 (s) ds 

1 + e~ 2T J-oo 

o T /2 poo 

= -^^(2e-T 2m e y2/{1+e ~ 2n e- s2 ^- 2 ^H 2m+1 (s)ds(4.37) 
1 + e~ /r Jy 

where to obtain the second equality the fact that the integrand is odd has been used. 

With x = (2iV) 1/2 + X/2 1 / 2 N 1 / 6 , making the change of variables s i-> (2N) 1 ' 2 + s and 
using the asymptotic expansion ( f4.17|) we see that 

x UN - 2m + l)\) l ' 2 e N2/H e Yt f°° e- st Ai(s) ds, u := 2m/iV 1/3 (4.38) 

JY+u 



For the x-dependent factor in ( |4.35| ), the asymptotic form of Jw(x)Cn(x) is given by 
( 4.18p with N — 1 replaced by N. Making use of the explicit value of h^, we thus find 



w e~ N f {tx - tY) Ai(X)e Yt J y e- st Ai(s)ds. (4.39) 



The asymptotic form of the quantity Si as defined by ( |2.28 ) is given by ( 4.21 ), since 
Si is the same for (3 = 1 and j3 = 4 initial conditions. Hence 

roc 

Sf\X,Y;t x ,t Y ) = / M(X + v)M(Y + v)e- v{tY ~ tx) dv (4.40) 

JO 

— Ai(X)e y * / e- st Ai(s) ds (4.41) 
2 Jy 
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Also, subtracting (p~3l"D divided by e'^'^PN 1 /* (and with t i-> t x - t Y ) gives 



Sf\X,Y-t x ,t Y ) = ~ [° Ai(X + v)M(Y + v)e~ v{tY - tx) dv 

--ki(X)e Yt / e~ st Ai(s) ds (4.42) 



Asymptotics of I(x,y;T x ,r y ) 

First note from ( [Hp , ( Q) , (p|) and (£J|) that 



w(x)e~^2^ k - 1 \^H 2k (x)e- 2T - (2A; - l)H 2k ^(x] 



.g~72fc 



-2- 2 ( fc - 1 )(|-( e - 2 / 2iy2fc _ 1 ( x )) _ Y!W^ (x)(e _ 2T _ 1} 



\dx^ y '> 2 

(c.f. ([4.271) ) • Use of ( [4.171) gives for the asymptotics 



§ N ->k if(2A) 1/2 + - ■ -L_] e iN-2 k ^t/N^rs 

_ e -i^/3t X7r l/4 2 7/4 2 -(iV-2fc-l)/2 (JV _ 2Jfe _ 1)11/2^1/12 ^_C eS tx Ai ( X + (4.43) 

s := 2k/ N. Similarly 



-$ 2fc _ 2 (y;r)e^- 2T = V«^)- 



r fc -i(0) vy ' 7 v vy/ 7T 1 / 2 (2fc - 1)! 

and so 



1 , / Y t 



-$7V-2 fc - 2 ((2A) 1 / 2 



o7iV-2fc-2 



r N /2-k-i(0) "- ZK - Z V ' 2 1 / 2 A 1 / 6 ' A 1 / 3 ' 

e -^/3t reaty7r -l/4 2 (JV-2fc-l)/2+l/4 (JV _ 2fc _ 1)1-1/2^-1/12 Ai( y + s y 

Multiplying together (|4.43| ) and ([4.44 ) and summing according to ( |2.21| ) we see that 

lf t (X,Y;t x ,t Y )= dse stY Ai(Y + s) — (e stx Ai(X + s))-(X^Y) (4.45) 

«/ do 

Asymptotics of D(x,y,T x ,T y ) 

Here we require the asymptotics of -R2m(-; r ) and i?2m+i( - ; t). The former is given by 
( 4.38|) . For the latter, we note from (f4.6|) that 



^ D /™. ^\„'72m + lT -"2m+l (y) 



r m (0) 2V7rr(2m + 2) 
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and thus according to ( [4.1 71) 



w(x) 



r/V/2-m(0) 



e Nr x e~ utx 



(4.46) 

7r l/4 2 (JV- 2m+ l)/2 + l/4 iv -l/12 Ai(x + u y 



z=(2Af) 1 /2 + x/2 1 /2 A ri/6 



2 v ^F(r(iV-2m + 2)) 1 /2 
Multiplying (|4.46|) and ( |4.38|) according to ( p.22| ) we see that 

£>(X, F; £ x , t Y ) = - e N / itx+tY) / e - uty Ai(y + u)e xtx / e - st Ai(s) - (X <-► y). 
4 Jo ix+B 

(4.48) 

5 Discussion 

5.1 Connection with (3 = 2 initial conditions 

We have calculated exact formulas for the dynamical distribution function (|2.8|) in the 
case of /3 — 1 and (3 = 4 initial conditions, with the final state corresponding to a (3 = 2 
equilibrium. If in these formulas we take the limit tx,t Y — > oo with — t Y fixed, then 
we would expect no memory of the initial state to be retained, but rather the distribution 
functions corresponding to a (3 = 2 initial condition (i.e. the equilibrium state) to result. 

Consider first the hard edge. Inspection of ( |3.31D , ( |3.37| ), ( |3.49| ) and ( |3.50D shows that 
in this limit 

S^ d (X, Y; t x , t Y ) ~ 5 4 hard (X, y, ; t x , t Y ) 

~ S 2 hard (X, y, ; t x , t Y ) := ~ jT 1 J Q (v^X) J a (VwF) e "«(**-*> , (5.1) 
S\ hard (X, Y; t x , t Y ) ~ S\ hard (X, y, ; t x , t y ) 

~ ^ iard (X, y ; t x , t Y ) := -- ^ J a (v^x) J a (v^l 7 ) e -«(*--^) d M , (5.2) 



while from ( ggg ), ( ggg ), ( gggp and (|337| ) we see that 



e (tx+ty) /h-d^ y. ^ ^ _ e (t,c+tr) y. ^ £y) ~ , (5.3) 

e -(tx+tr) 7jhard( X? y. ^ ^) _ g-(tx+*r) Tjhard^ y. ^ ty) ~ . (5.4) 



Substituting these values in Q2.24 ) and then using ( |2.26 ), and the formula 



Pf(ZX) = (det(ZX)) 



1/2 



we see by interchanges of particular rows and columns that ZX is equivalent to a block 
matrix of the form 

(On+m A \ 
~A On+m J 
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where 



.4 



Chard / vi. 1 ) vi 1 ) . + + 
D 2 l A i > A fc 



Chard / Y"M y( 2 ).+ + 

^2 ^ A j 5 A fc 5^15^2 

Chard 

/ v (2) v (2) ,. . 
^2 y^j > A fc ) Z; 2,t2 



Chard f V (2) v (l) , . 

^2 l/S' > A fc ) J 2,r 
Thus the limiting form of the dynamical distribution is given by 

nt s (r {1) r (l). T . r ( 2 ) T (2). _ \ _ A pf A 

P(n+m) \X 1 ; • • • ) x n ) r l> x l ) • • • ) x m i T 2) — Uei rv . 



(5.5) 



(5.6) 



The analysis in the soft edge case is very similar. The formulas ( |5.3|) and (5.4) again 
apply, as do the formulas ( |5TT| ) and ( |5.2| ) with S)^ xA , S\ wA replaced by 



S s 2 oit (X,Y;t x ,t Y ) := / M(X + v)M(Y + v)e~ v{tY - tx) dv , 

Jo 

Sf\X,Y-t x ,t Y ) := - f M(X + v)ki(Y + v)e- v{tY ~ tx) dv. 

J — oo 

We remark that in the case n — m — 1, the resulting formula for p( n+m ) implied by ( |5.6| ) 
and (|5.5| ) in both the hard and soft edge theories agrees with that computed by Macedo 



(5.7) 
(5.8) 



3 0. 



5.2 Connection between the soft and hard edges 

For the hard edge distribution function of the Laguerre ensemble at (3 = 2 and r = 0, it 
is known that by taking the limit a —>■ oo and introducing new scaled coordinates, the 
distribution function of the soft edge of the Gaussian ensemble results . Specifically, 
the new scaled variable x is related to the scaled variable X of the hard edge by 



X 



2a(a/2) 1/3 x 



(5.9) 



and the distribution functions are related by 

Jim (-2a(a/2) 1 / 3 )"p^ d (a 2 - 2a(a/2) 1 ' 3 x 1 , 

= P(n) 



a 2 - 2a(a/2) 1/3 x n ) 



(5.10) 



As noted in [IB], |T7| this property should persist for general (3. Further, with appropriate 
scaling of r, the result would be expected to generalize to the dynamical distributions. 

Let us first check this latter point for the dynamical distribution with a (3 = 2 initial 



condition ( |5.6|) . For this purpose we require the asymptotic expansion |L8] 

Ja{x) r 



2\ 1/3 A .(2^ /3 (a-xY 



X 



X 



1/3 



(5.11) 



valid for a and x large such that the argument of the Airy function is of order unity. 
Changing variables u = 1 — w(2/a) 2 ^ 3 in the integrals of (|5.1j) and (|5.2j ) and introducing 
the scaled coordinates ( |5.9| ), application of ( |5.11| ) shows 



J a (v^) 



2 v 1/3 

- ) Ai(x + w ) 
a J 



(5.12) 
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Thus if we also scale t according to 

t x >-> {a/2f/H x 

then we see that 

-2a(a/2) 1 / 3 ) S 2 hard (a 2 - 2a(a/2) 1 ' 3 x,a 2 - 2a(a/2f' 3 y- (a/2) 2 ^t x , (a/2) 2 'H 

-(a/2) 2 / 3 (t x -t y )nsolt( f f \ 

u 2 v x ' vi b xi y) i 



(5.13) 



~ e 



(-2a(a/2) 1/3 ) S 2 hard (a 2 - 2a(a/2) 1/3 x, a 2 - 2a(a/2) 1/3 y; (a/2) 2 / 3 **, (a/2) 2/ \ 

Since, as seen from ( |2.24| ), the factors e~^ a//2 ^ 2/3< - t:c_t! '- ) do not affect the value of p( n+m ), we 
thus have the dynamical extension of ( [5.10 ), 

Jim (-2a(a/2) 1/3 ) n+m p^ d m) (a 2 - 2a(a/2) 1/3 x 1} . . . , a 2 - 2a(a/2) 1/3 x n ; {a/2) 2/3 t x ; 

a 2 - 2a(a/2) 1 / 3 y 1 , . . . , a 2 - 2a(a/2f/ 3 y m - (a/2) 2 /%) 

= P\n+m)( x li ■ • • i X n] tx] V\i ■ ■ ■ > Vm\ ty) (5-14) 



Consider next (3 = 1 initial conditions. Inspection of ( |3.31| ), ( |3.37|) , (|3.33|) and (|3.38 



shows that the analysis used above for (3 = 2 initial conditions again suffices to establish 
the connection formula. For example, consider the formula ( |3.33| ) for Jj iard . For large a 
we see that 

I\"\X,Y;t x ,t Y ) 



-uti 



U 



-J a (VuX) du-(X <-> Y). (5.15) 



Use of the asymptotic expansion (|5.12 ) ( after making the appropriate change of variables 
in the integrals) shows 

-2a(a/2) 1 / 3 ) I*"* (a 2 - 2a(a/2) 1 / 3 x, a 2 - 2a(a/2) 1 / 3 y; (a/2) 2 / 3 ^, (a/2 f 3 t y 



~ e 



Similar formulas hold connecting S'j iard , S'j iard and D\™ A to S s ° lt , Sf h and Df lt respec- 
tively, thus showing the formula ( |5.14| ) holds for (3 = 1 initial conditions. Repeating the 
analysis for (3 = 4 initial conditions leads to the same conclusion. 



5.3 Connection with the bulk 

It is known |HJ that after appropriate scaling of the variables, the static distributions at 
the hard and soft edges tend to the bulk distributions at large distances from the edge of 
the system. The required scaling is 

X ^ a + 2ix\fci px , X \— > — (a + irpx/\/~a} , (5.16) 
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with a — > oo, for the hard and soft edges respectively. Here we will show that this feature 
persists for the dynamical distributions. 

Consider first the hard edge with (3 — 2 initial conditions. Changing variables u i— > u 2 
and making use of the asymptotic expansion 



x-^oo 

in the formulas (15.11) and (15.2), we see that 



( 2 \ 1/2 
— ) cos(x — na/2 — 7r/4) 



(5.17) 



lim 2icy/a p S , o ard (a + 2'n^fa px, a + 2-Ky/a py; t x , t v ) 

a — >oo \ / 

~ p cos(nup(x - y)) e ~ u2{tx ~ ty) du , 
Jo 

lim 2'Ky/a p Sj ard [ct + 2-K^fa px, a + 2iry/a py; t x , t v 



COS 



(-nup(x-y))e u2{tx ty) du. 



(5.18) 



(5.19) 



Substituting these formulas in ( |5.5| ) and ( |5.6| ) and also making the change of scale t i— > 
(7Tp) 2 t/2 reclaims the known formula for the bulk dynamical distributions 

In fact the asymptotic behaviour (|5.18 ) and (|5.19| ) of S^ ard and S^ ard for (3 = 2 initial 
conditions persists for (3 = 1 and (3 = 4 initial conditions. This follows from application 
of ( |5.27| ) in the exact results Q3.31J) , ( |3.37| ) and ( |3.49| ), ( |3.50| ). It remains to consider Ip 
and Dp. For (3 = 1 initial conditions application of ( 5.17 ) in (|3.33 ) and ( 3.3S ) shows 

lirn^ 2np I^ rd (a + 27ry / a px, a + 2it\fa py; t x , t^j 

roc Q—u 2 (t x +t y ) 

= pi du sm(irup(y — x)) , (5.20) 

Ji u 

lim 2npa D\ ard (a + 2-n^fa px, a + 2ix\fa py; t x , t y j 

= p[ duue u2{tx+ty) sm(irup(y - xj) , (5.21) 

JO 

while for (3 = 4 initial conditions application of ( 5 . 1 7|) in ( |3.54| ) and ( 3.57|) shows 
lim 2npa ii ard [a + 2ny/a px, a + 2iry/a py; t x , t v ) 

du u e~ u (tx+ty) sm(-Kup(y - x)) , (5.22) 
lim 27rp.D$ ard (a + 2ir\/a px, a + 2-Ky/a py; t x , t v 



p I du 



U 2 (t X +ty) 



u 



sm(irup(y — x)) . 



(5.23) 



After scaling t \— > (iip) 2 t/2, the known formulas for y; t x , t y ) with (3 = 1 and 

(3 = 4 initial conditions J2t], ^TJ are reproduced by these results. 

In the case of the soft edge connecting to the bulk, similar formulas are found to hold. 
These formulas are derived from the asymptotic expansion 



Ai(-x) 



TF> — 777 COs(2x 3/ 73 - 7t/4) . 

x^oo ^1/2^,1/4 v ' I ' 



(5.24) 
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Consider for example J S'| oft as specified by (|5.7| ). We first change variables v i— > a(l — u) 
for < v < a, and replace the Airy functions in the integrand according to (|5.24j) . Noting 
that 

— I (a + ixpx/ \fa — t>) 3 / 2 — (a + npy/y/a — v ) 3 ^ 2 } ~ ■Kpu 1 ^ 2 {x — y) , 

we see that if we also scale t according to 

th^t/a (5.25) 

then 5*2 has the asymptotic form 

lira (— irp/y/a) (—(a + npx/y/a), —(a + Tcpy/y/a);t x /a,t y /a) 

~ e~ {tx ~ ty) p / cos(tcup(x - y)) e~ u2{tx ~ ty) du . (5.26) 
J o 

Apart from the factor of e~^ tx ~ tv ^ this is the same result as (|5.18|) . A similar calculation 



shows the asymptotic form of Sf * s given by ( p.iyp except for an additional factor 
of e~( tx ~ ty \ This shows that this additional factor cancels from the formula (|2.24 ) for 



the distribution functions. The final result is that the bulk dynamical distributions are 
reclaimed from the soft edge results for each of the three distinct initial conditions. (As 
in the hard edge case a final rescaling t i— > (np) 2 t/2 is required). 



5.4 Asymptotic form of pT 1+1 \(X,Y;t) 

The dynamical density-density correlation function between an eigenvalue with value X 
initially, and an eigenvalue with value Y and parameter r is defined in terms of the 
distribution function (|2.8|) by 

pfi+i) (x, y; r) = p (1+1) (x; 0; y; r) - p (1) (x; 0) p {1) (y; r) . (5.27) 

According to (2.24) and fl2.26| ) we have that in the finite system 



pj l+1) (x,y;r) = - [S(x, y; 0, r) S(y, x; r, 0) - I(x, y; 0, r) D(x, y; 0, r)) . (5.28) 
For the Laguerre ensemble at the hard edge we define the corresponding scaled correlation 

pfi+T) d by 

T hard / v i/ - ,\ i- / 1 \ <j< / ^ ^ 

(A, r ; t) = jim ^— j p (1+1) ^— , — ; — ^ 

^J ard (X, F; 0, t) 5j ard (F, X; t, 0) - iJ ard (X, Y; 0, t) Z^ ard (X, F; 0, tj) ,(5.29) 



where /3 = 1, 2 or 4 depending on the initial conditions, while for the Gaussian ensemble 
at the soft edge the scaled correlation P(i+i) is defined by 

pl^(X,Y;t) = ^(2^^)^ 

= - (^ oft (X, Y; 0, t) Sf\Y, X; t, 0) - lf\X, Y\ 0, t) Df\X, Y\ 0, t)) . 

(5.30) 
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The asymptotic form of Y; t) for large \X — Y\ and t is of interest in applications 

of the parameter dependent Laguerre ensemble to transport in mesoscopic systems |15| . 
To deduce this asymptotic form, a hydrodynamical approximation to the corresponding 
Fokker-Planck equation was made which implied 



j o j( (4XY) p£$f (x 2 , Y 2 ; 1) cos(kX) cos(kY) dX dY ~ q U e ~^\ ( 5 .31) 



(eq. (80) of |15|] with irpdu 2 /^ = 2t). On the other hand we can easily check that 



1 1 



Re I t + t -77 I = -2 / dk ke~ 2kt cos(kv) cos(ku) . (5.32) 

1 [y + u + 2it) 2 (v-u- 2it) 2 1 ' 



OO 



Taking the inverse transform of this result and comparing with ( |5.31| ) shows that the 
latter formula is equivalent to 

(«»f l5 ? {X\ y> ; t) --L R e ( {X + Y \ 2tt)2 + (x J +2lt) ) (5-33) 

for the leading non-oscillatory behaviour. 

The prediction (|5.1| ) can readily be checked on the exact formula ( |5.29|) with (3 = 1, 2 
and 4 initial conditions. The simplest case is (3 = 2 initial conditions, when the product 
jhard£)hard j n ( 5 29|) vanishes. To analyze S < 2 ard and S"2 ard we make use of the asymptotic 
expansion ( |5.17 ), and furthermore expand the integrands about their endpoint u = 1. 
Integration by parts then gives 

t ( J{X-Y) i((X+Y)-n a /2-7T/2)\ 

So ard fx 2 , Y 2 ] 0, t) ~ Re I— I + f_ , 

2 V ' ' '7 27r(Xy)V2 \-i(X -Y)-2t —i(X + Y) — 2t J 

(5.34) 

e -t / gi(X-r) e i((X+y)-7ra/2-7r/2) ^ 



^ d (F 2 ,X 2 ;t,0) ~ 27r(XF) i/2 Re ( v , (X _ y) _ 2t + ,( X + F)-2t 



(5.35) 



and these formulas when substituted in ( 5.29Q imply the result ( 5.33| ) with (3 = 2. In the 
case of (3 = 1 and (3 = 4 initial conditions, similar analysis shows 

jhard ^2 ; y2. ^ ^hard ^ y 2. ^ 

~ XpS^ d (X,Y;0,t)S^ d (Y,X;t,0) 

~ X "IXY 2^ Re ( (X + Y + 2U) 2 + (X - Y + ' ^ 

where xp = ~ 1 f° r = 1 an d X/3 = 1/2 for /3 = 4, and the asymptotics refer to the 
leading non-oscillatory term. Thus ( |5.33| ) again remains valid. 

Let us now turn our attention to the same asymptotic limit for the soft edge. There 
it is known that the asymptotic form ( |5.33| ) for X, Y — > oo is valid in the static theory 
(t = 0) [|22| . It turns out that the exact results for ( |5.30|) satisfy (|5.33j ) for X, Y — ► oo 
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and t — > oo, as just demonstrated for the hard edge correlations. Thus, for (3 = 2 initial 
conditions, use of the asymptotic expansion ( |5.24j ) in ( |5.7| ) and (|5.8|) shows 



S so{t ( X 2 ,-Y 2 -Q,t 



j / e 2i(X :i +Y 3 )/3-iv(X+Y)-m/2 ^{X^-Y^/^-iviX-Y)' 

Re — — — ■ + 



2tt(XF) 1 /2 ^ i(x + Y)+t i(X-Y)+t 



/ e 2i(X 3 +Y s )/3-iv(X+Y)-wi/2 e 2i(X 3 -Y 3 )/3-iv(X-Y) ' 

Re — — — ■ + 



2tt(XF) 1 /2 ^ i(X + Y)-t i(X-Y)-t J' 

and these formulas when substituted in ( |5.3C| ) (with If Df = 0) gives the asymptotic 
form ( |5.33| ) (although t therein must be rescaled by replacing 2t with t). The same 
conclusion is reached for (3 = 1 and 4 initial conditions, where we find the analogue of the 
formula ( |5.36| ). 

As pointed out in fll5| , a consequence of the asymptotic form (|5.33|) is a universal 



formula for the time displaced covariance of two linear statistics 

N N 
3=1 J=l 

Thus, if in the general formula 

poo poo 

Cov(A , B t ) = dxa(x 2 ) dyb(y 2 )(Axy)pf 1+1) (x 2 ,y 2 ;t) 

we suppose a(x 2 ) and b(x 2 ) are made slowly varying by writing 
a(x 2 ) h-> a((x/a) 2 ) , b(x 2 ) i-> b((x/a) 2 ) , 
where a>l, and if we scale t according to 

1 1 — > oc t , 

then a simple change of variables and use of ( |5.33| ) shows that for a — ► oo 



Cov(A ,B t ) — J^ dxa(x 2 )J Q dy b(y 2 )Re f — - — 



+ 



+ 2it) 2 (X-r + 2^)^ 
dkkd{k)~b{k)e- 2kt , (5.37) 



7T 2 /3 Jo 

where the equality follows from ( |5.32| ), and 

("OO 

a(fc) := / a(x 2 ) cos(/cx) eke 
and similarly the meaning of 6(fc). 
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5.5 Empirical density of eigenvalues at the soft edge 



The dynamical distribution functions at the soft edge calculated in Section 4 can be real- 
ized by parameter dependent random matrices of the type (|1.3j) . Thus we can numerically 
construct random matrices and empirically compute eigenvalue distributions for compar- 
ison against the theoretical prediction. The simplest distribution to compute empirically 
is the density (one-point function). Here we will undertake such a calculation for the 
orthogonal to unitary transition at the soft edge. 

As deduced from (|2.24|) , (|2.26|) and (|3.31|) the theoretical prediction for the scaled 



density at the soft edge of the orthogonal to unitary transition is 

f$(X;t) = Sr ft (X,X;M) 

= / (M(X + v)) dv + - / e~ stx Ai(X - s)ds. (5.38) 
Jo v ' 2 Jo 



The theory of Section 1 tells us that the parameter dependent random matrices ( |1.3|) , for 
an appropriate X^ ' have, in the infinite dimension limit, a scaled density at the soft edge 
given by ( |5.38| ). Since we are considering the orthogonal to unitary transition, must 
be chosen as a real symmetric matrix with joint p.d.f. for the elements 

A / 1 \ 1 / 2 y( ' 2 /9t-t 1 x^ 2 

^ (o) ) = n(^) e~ x - 1 . 

j=l j<k 

The joint distribution of the elements of X, P(X; r) say, is obtained by integrating over 
xf) and X$. Thus 

A' 



/OO P oo 

dX® n / dX^P(X^)P(X^;X;r) 



j<k 

N -X 2 



Yl 6 n /( 1 + e -Q 2 r -2(ReX lk ) 2 /a+e- 2 -) r ~2amX lk ) 2 /(l-e- 2 -) 



j=i x /n(l + e- 2r ) j<k n vl 



e 



-At 



(5.39) 



which tells us that the diagonal and the real and imaginary parts of the upper triangular 
elements are all independent, having Gaussian distributions with mean zero and variances 
2/(1 + e~ 2r ), 4/(1 + e~ 2r ) and 4/(1 — e~ 2r ) respectively. Hence, for a given value of iV 
and r such matrices are simple to generate numerically. For each such random matrix the 
eigenvalues Xj are computed and scaled according to A i— ► (Xj — (2A r ) 1 / 2 )2 1 / 2 N 1 ^ 6 . For 
large enough and with r related to t by ( |4.9| ) the corresponding empirical density should 
approach ( |5.38D . The situation for a relative small value of N (recall the occurence of 



iV 1//6 ) is illustrated in Figure 1, where good agreement with the infinite dimension limit 
is observed. 

Also of interest is the possibility of providing a realization of the hard edge results. 
Unfortunately here it seems that the initial conditions used do not precisely correspond to 
a situation in which X in (|1.6|) can be generated from a formula analogous to ( |5.39| ) . The 



problem is that if X^ is a real n x m matrix with independent Gaussian entries, then 
the square of the eigenvalues of X^X^ has p.d.f. ( |1.2| ) with a = n — m — 1. However in 
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density 




Figure 1: Empirical histogram of the scaled soft edge density for 500 parameter dependent 
random matrices of the form (|5.39|) with r = .05 and N = 100, compared against the 
infinite dimension form ( |5.38| ). 



deducing ( |1.4|) from ( |1.6D in the case of transitions to unitary symmetry ((3 = 2) we must 
take a = n — m, which is thus incompatible with the sought initial condition. 



5.6 Distribution functions in the initial state 

With m = and T\ — > the dynamical distribution ( |1.8| ) reduces to the n-point distribu- 
tion function in the initial state. At (3 = 1 the formulas (gjj) , ( |2T30D , ( [2"3lD and (^33|) 
reclaim the well known [JTSJ, |TD|, || ^3] general formula 



P(n){xi, ...,x n ) = Tdet 



Si[Xj, Xk) 

Q x . {.Xj 5 Xk 



§ / sgn(x k - y)Sx{xj, y) dy - ±sgn(x fe 



(5.40) 

where Si(xj, Xk) := 5*(xj, x^; 0, 0) as specified by ( |2.20p with {Rj} corresponding to the 
(3 = 1 initial condition. We remark that the integral in ( |5.40|) can be rewritten according 
to 



- J sgn(x fc - y)S 1 (x j , y) dy 



Si(xj,y) dy. 



At (3 = 4 these same formulas (with ( |2.33| ) replaced by ( |2.34| )) give that 
p( n )(xi, ...,x n ) = Tdet 



S^(xj,Xk) q x ^ S\ {xki Xj / 



Sxl S i{Vi x h) dy S 4 (x k} Xj) 



(5.41) 



where S^x.y) := S(x, y; 0, 0) as specified by fl2.20|) with {-Rj} therein corresponding to 
the (3 = 4 initial condition, again in accordance with the known result. In deriving ( |5.41| ) 
from (|2.24j ) we assumed that none of the points coincident; this is necessary 

because of the delta functions in ( |2.7| ) with (3 = 4. 

We can take the limit tx,ty — > in the formulas obtained for Si se (X, Y; tx, ty) 
(here edge denotes hard or soft) obtained in Section 3 and 4 to obtain the explicit form 
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of SI ge (X, Y) in the scaled limit at the hard and soft edges. Thus from (|3.31|) , ( |3.49|) 
( ggg ) and dO0| ) we find 



where 



= ^y^^j// J a+1 («)d W -l) (5.42) 

4v^ "'o 

S 4 hard (X,y) = K^\X,Y)-^^J a ^i{X 1 ' 2 ) J a+ i(s)ds (5.43) 

1 . POO . 

Sf^X.Y) = X soft (X,Y) + -Ai(Y)(l - / Ai(t)dt) (5.44) 

I /'OO 

^ oft (X,y) = X soft (X,Y) - -Ai(X)J Ai{s)ds (5.45) 
XWj^X^JjyW) - Y^J a+ i(Y^)J a (X^) 



K hara (X,Y) :-- 



2(X - Y) 



«- M := fflMWl, (5.46) 

X — Y 

In obtaining these formulas we have used fl5.2| ) and ( |5.7| ) together with the known facts 
that 

sJ? ard (x, y ) = K h&rd (x, y ) , ^ oft (x, y ) = x soft (x, y ) 



to rewrite the first of the integrals in the expressions (|3.31| ), ( |3.49|) , (|4.22|) and ( |4.40|) 



Also, Bessel function identities have been used in obtaining the second term in ( |5.42|) and 
( |5.43| ), while in the second term of ( |5.44|) use has been made of the definite integral 



oo 

Ai(x) dx = 1. 

-oo 



We remark that in previous studies || [24] expression have been obtained for (|5.42 ) and 



(|5.43|) which are of a more complicated structure. Also the expressions obtained in [16| 



for ( |5.44j ) and ( |5.45| ) are now seen to be in error. 



Some caution must be exercised in applying these results at f3 — 4. In particular, 
the one body weight function used in (|2.7|) at (3 = 4 is that corresponding to (3 = 2 in 



(|1 . 1| ) and (|1.2| ). To obtain from the above results expressions corresponding to the weight 
functions as written in ( |1 . 1| ) and ( |1.2| ) we must make the replacements 

5 4 hard (X, Y) >-> 2S 4 hard (4X, 4Y) , Sf ft (X, Y) i-> -^S? ft (2 2 / 3 X, 2 2 / 3 Y). 



a _2a' 4 v">* y 2 i/ 3 ' 
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Appendix A 

In the = 0, the formula ( |3~9|) for 52 at (3 — 1 reads 

S 2 (x,y;0,0) = v^( ^4ff ~ ^)^-^) (C N „i(y)-(N-l)R N - 2 (y)). (A.l) 
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Use of the explicit formulas flO| ) and (|3.2|) gives 



C N -M -(N- l)R N - 2 {y) = -(N - l)\(L a N _M - ^Sr_i(v)) = (N - l)^L a N (y), 



dy 



(A.2) 



where the second equality follows from an appropriate Laguerre polynomial identity 

For the x-dependent factor in (|A.1|) , we require an explicit formula for $tv-2(^;0) 
Now ( |2.31| ) substituted in fl2.16| ) gives 



®k(x;0) = / sgn(x - y)^w(y)R k {y). 



Making use of the explicit formulas ( |3.2| ) and ( |3.3| ) then gives 



— — / sgn(x - y)y a/2 e-y/^L%_ 1 (y) dy. 
i + N ) Jo dy 



^2-1(0) AT(a + N)Jo ~°" v a,a dy 
Integrating by parts, and making use of the Laguerre polynomial identity 

(y - a)L a N ^(y) = y^-L^y) - N(L%{y) - 2&-i(v)). 



we see from this that 



r N /2-i(0) 



W(X) 



dy 
-C N Ux 



h 



N-l 



N 



4T(a + N) 

Substituting flA~2] ) and in flO]) shows 

TV! 



j o sgn(x - y)y^ 2 - 1 e-y/ 2 (L a N (y) - L%_ x 



dy. 



(A.3) 



S 2 (x, y; 0, 0) 



AT{N + a) 



w(y) 



x(-L%(y)) / sgn(x -y)y a / 2 - 1 e^ 2 (L%(y)-L a N _ 1 (y))dy, 



(A.4) 



-dy " ' " ' Jo 

which is the formula recently given by Widom ||. 

Consider next the expression ( |3.42j ) for S 2 at (3 = 4 with r x = r y = 0. To simplify this 
expression, we note that the formula in ( |3.39|) for R 2 j(x) can be written in the simplified 
form 



R 2j (x) 



a + 2j 



±—(Cv +1 (x) + \{2j + l)\x^ffi £t^e^(±IwMt)) *)-(A.5) 



This formula can be checked by multiplying both sides by x a ^ 2 e~ x ^ 2 and differentiating. 
(We also verify that the x — > 00 behaviour of both sides agrees.) Now, according to the 
definition ( |2.16| ) of and the second formula in ( |2.32j ) 



d 



$jv- 2 (x;0) _ 2 
rjv/ 2 -i(0) r N/2 -i(0)dx 



(x a ' 2 e~ x ' 2 R N . 2 i 



x 
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Substituting ( A.5 ), and making use of the differentiation formulas 



y^K(y) = nL a n (y) - (n + ajL^y) = (n + l)L a n+1 {y) - (n + a + 1 - y)L»(y) 



and the fact that 

we readily find 

$jv_ 2 (x;0) 
rjv/2-i(0) 



(a + N- l)r N/2+1 (0) = h N -! = T(N)T(N + a) 



w(x) 



C N -i{x) 
h N _i 



|r|r| J (L M x)-^_ lM ) 



T(AT + a) 



which thus simplifies the x-dependent factor in ( 3.42 ). For the ^/-dependent factor, we see 
from ( TOp that 



V^(CWy)-(a + iV-l)i?Ar- 2 (y)) =-\{N -\)\ £r*e-W±L%{t)dt. (A.7) 
Multiplying ( A.6|) and ( |A.7 ) shows 



if (A* + a) x Jo at 



x Jo 

in agreement with the formula given recently by Widom 



(A.8) 
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